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I. INTRODUCTION 

Theories unifying gravity with other interactions sug- 
gest a possibility of temporal and spatial variation of 
the fundamental constants of nature (see e.g. review 
0| where the theoretical models and results of measure- 
ments are presented). There are hints for the variation 
of the fundamental constants in Big Bang nucleosynthe- 
sis, quasar absorption spectra and Oklo natural nuclear 
reactor data. However, a majority of publications report 
only limits on possible variations. For example, compari- 
son of different atomic clocks gives limits on present time 
variation of the fundamental constants. 

A large fraction of the publications discuss the vari- 
ation of the fine structure constant a = e 2 /he. The 
hypothetical unification of all interactions implies that 
a variation in a should be accompanied by a variation 
of the strong interaction strength and the fundamental 
masses. For example, the grand unification model dis- 
cussed in Ref. |2| predicts the quantum chromodynamics 
(QCD) scale Aqcd (defined as the position of the Lan- 
dau pole in the logarithm for the running strong coupling 
constant) is modified as SAqcd/Aqcd ~ 34 8a/ a. The 
variation of quark and electron masses in this model is 
given by 8m/m ~ 70 5a/a, giving an estimate of the 
variation for the dimensionless ratio 

S(m q /A Q cp) _ 35 8a 
(m q /A QC D) a 

The coefficient here is model dependent but large val- 
ues are generic for grand unification models in which 
modifications come from high energy scales; they appear 
because the running strong-coupling constant and Higgs 
constants (related to mass) run faster than a. If these 
models are correct, the variation in quark masses and the 
strong interaction scale may be easier to detect than a 
variation in a. 

One can only measure the variation of dimensionless 
quantities. We want to extract from the measurements 
the variation of the dimensionless ratio m q / Aqcd where 
m q is the quark mass (with the dependence on the renor- 
malization point removed). A number of limits on the 
variation of m q / Aqcd have been obtained recently from 
consideration of Big Band nucleosynthesis, quasar ab- 



sorption spectra and the Oklo natural nuclear reactor, 
which was active about 1.8 billion years ago 

Karshenboim G] has pointed out that measurements of 
ratios of hyperfine structure intervals in different atoms 
are sensitive to variations in nuclear magnetic moments. 
However, the magnetic moments are not the fundamental 
parameters and can not be directly compared with any 
theory of the variations. Atomic and nuclear calculations 
are needed for the interpretation of the measurements. 
Below, we calculate the dependence of nuclear magnetic 
moments on m q / Aqcd by building on recent work and 
incorporating the effect of the spin-spin interaction be- 
tween nucleons. We obtain limits on the variation of 
m q l Aqcd from recent experiments that have measured 
the time dependence of the ratios of the hyperfine struc- 
ture intervals of 133 Cs and 87 Rb 0, 133 Cs and *H 0, 
m Yb+ and 133 Cs 0, 199 Hg+ and X H [nj, the ratio 
of the optical frequency in 199 Hg + to the hyperfine fre- 
quency of 133 Cs [l2j, the ratio of the optical frequency 
in *H to the hyperfine frequency of 133 Cs ^3|, and the 
ratio of the optical frequency in 171 Yb + to the hyperfine 
frequency of 133 Cs ^4j. It nas been suggested in Ref.fi1| 
that the effects of the fundamental constants variation 
may be enhanced 2-3 orders of magnitude in diatomic 
molecules like LaS, LaO, LuS, LuO. Therefore, we also 
present the results for 139 La. 

During the calculations, we shall assume (for no- 
tational convenience) that the strong interaction scale 
Aqcd does not vary and so we shall speak about the 
variation of masses (this means that we measure masses 
in units of Aqcd)- We shall restore the explicit appear- 
ance of Aqcd m the final answers. 

The hyperfine structure constant can be presented as 

A = const x (^) [a 2 F rel (Za)] ( M — ) (2) 
nr m p 

The factor in the first set of brackets is an atomic unit 
of energy. The second "electromagnetic" set of brackets 
determines the dependence on a and includes the rela- 
tivistic correction factor (Casimir factor) F re i . The last 
set of brackets contains the dimensionless nuclear mag- 
netic moment /i (that is, the nuclear magnetic moment 
M = fi[eh/2m p c] ) and the electron and proton masses 
m e and m p . We might also have included a small correc- 



2 



tion due to the finite nuclear size but its contribution is 
insignificant. 

The ratio of two hyperfine structure constants for dif- 
ferent atoms will cancel out some factors such as atomic 
unit of energy and m e /m p , and any time dependence 
falls on two values: the ratio of the factors F re i (which 
depends on a) and the ratio of the nuclear magnetic mo- 
ments (which depends on m q / Aqc d) ■ 

For the F re i component, variation in a leads to the 
following variation of F re i [TjJ 



5F r , 



Sa 



F 



rel 



(3) 



and one can use the s— wave electron approximation for 
F re i to get 



K r 



(Za) 2 (12 7 2 - 1) 
7 2(4 7 2 _ 1} 



(4) 



where 7 = y 1 — (Za) . However, numerical many- 
body calculations [16j give more accurate results, with 
a slightly higher value of K re i than that given by this 
formula. A comparison is shown in Table [I] 

The other component contributing to the ratio of two 
hyperfine structure constants is the nuclear magnetic mo- 
ment /i. Theoretical values for /i in the valence shell 
model are based on the unpaired valence nucleon and are 
given by Schmidt values 



for j = I 



2(i+ 



jy \-g s + (2j + 3)g{\ for j = I - § 



(5) 



The orbital gyromagnetic factors are gi = 1 for a va- 
lence proton and gi = for a valence neutron. The spin 
gyromagnetic factors are g s (= g p ) = 5.586 for protons 
and g s (= g n ) = —3.826 for neutrons. These g-factors 
depend on rn„ / Aon n and previous work exploring this 
dependence [171 II81 ] is summarized below. We then use 
these results to consider the more realistic situation of \x 
having both a valence nucleon contribution and a non- 
valence nucleon contribution due to the spin-spin inter- 
action. 



II. VARIATION OF MAGNETIC MOMENT 
WITH VARIATION IN QUARK MASS 

A. Variation in fi using valence model magnetic 
moment 

As a preliminary to our results and as a comparison for 
evaluating the effects of our calculations, we include the 
results of work previously done in this area 

[HG1. This 

work was essential to our results as the authors calcu- 
lated the variation in the neutron and proton magnetic 
moments {[i n and fi p ) with the variation in m q j Aqcd 
using chiral perturbation theory. 



As mentioned above, the g-factors depend on 
m q / K.QCD- The light quark mass m q = {m u + m,d)/2 « 
5MeV and in the chiral limit m u = m^ = 0, the nu- 
cleon magnetic moment remains finite. Thus one might 
assume that corrections to the spin ^-factors g p and g n 
are small. However, the quark mass contribution is en- 
hanced by 7r-meson loop corrections to the nuclear mag- 
netic moments, which are proportional to 7r-meson mass 
771^ ~ \fm q KQCD- Since m^ = 140 MeV, the contribu- 
tion can be significant. 

Full details of these calculations are given in Ref 
[T7L ITsj . They give the following results, which relate 
variations in fi n and fi p with variations in light and 
strange quark masses (m q and m s ): 



MP 


0.087 5m « 

m q 


(6) 


MP 


8m* 

0.013 

m s 


(7) 




0.118 Sm « 

m q 


(8) 


<W 

Mr. 


<5m a 

+ 0.0013 s 

m s 


(9) 



Using these relations and the valence model approxi- 
mations for fi, we can obtain expressions of the form 



6 ft, 



5m„ 



fx Tn q i,fj s 

Hence for nuclei with even Z and a valence neutron 
Sji _ Sg„ 
M g n 



(10) 



-0.118 q - + 0.0013 



For valence protons, the orbital gyromagnetic factor gi 
also has an effect. Thus for 133 Cs with its valence proton 
and j = I — I , 

^ = 0.110^+0.016^ 



while for 87 Rb with its valence proton but j = I + | , 

^ = _ .64^- 0.010^ 
fi m q m s 

These results can be presented using the ratio of the 
hyperfine constant A to the atomic unit of energy E = 
m e e /ti 2 by defining the parameter V through the rela- 
tion 



5V_ _ 5{A/E) 
V ~ A/E ' 



(11) 



The values for n q and n s in the results for 5/i/fj, can 
then be combined with the corresponding values of K re i 
in Table to give results of the form: 



V(M) = a 



A 



QCD 



)K 3 
— (12) 
m p 
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TABLE I: Variational factor K re i and k, values for various atoms obtained using simple valence shell model (method A) as used 
in equation 11211 . The first row is given by equation The second row presents the results of the more accurate many-body 
calculations (see ref. (l6jL The numerical results marked by * are obtained by an extrapolation from other atoms. 



Atom 


X H 


2 H 


3 He 


87 Rb 


ln Cd 


129 Xe 


133 Cs 


171 Yb 


"»Hg 


K re i (analytical) 











0.29 


0.53 


0.71 


0.74 


1.42 


2.18 


K re i (numerical) 











0.34 


0.6* 


0.8* 


0.83 


1.5* 


2.28 


Kq 


-0.087 


-0.020 


-0.118 


-0.064 


-0.118 


-0.118 


0.110 


-0.118 


-0.118 




-0.013 


-0.044 


0.0013 


-0.010 


0.0013 


0.0013 


0.016 


0.0013 


0.0013 



The calculated values of K re i and k to be used in the 
expression for V(M) for various atoms are summarized in 
Table [I] The factor m e /m p will cancel out when a ratio 
of hyperfine transitions is used. It will, however, survive 
in a comparison with optical and molecular transitions. 



B. Variation in ^ incorporating the effect of 
non-valence nucleons by using the experimental 
magnetic moment 

The results of the previous section were used to cal- 
culate the variation of /i with m q / Aqcd based on the 
single particle approximation for fi (one valence nucleon) 
within the shell model. That is, it was assumed that 
the dimensionless nuclear magnetic moment is given by 
/i = g. p (s z ) +(l z ) for a valence proton, and /i = g n (s z )° 
for a valence neutron. Here, g p and g n are the spin gy- 
romagnetic factors for free protons and neutrons respec- 
tively, (l z )° = jz — (s z )°, and (s z )° is the spin expectation 
value of the single valence nucleon in shell model: 



(O = 



2(J + 1) 



for j = I + | 
for j = I - | 



(13) 



However, it is well known that this theoretical value 
is only an estimate of \x and the magnetic moment of 
the valence nucleon tends to be offset by a contribution 
from the core nucleons. An empirical rule is that the spin 
contribution of a valence nucleon should be reduced by 
40% to obtain a reasonable value for the nuclear mag- 
netic moment. This reduction may be explained by the 
contribution of core nucleons, which should be negative 
since proton and neutron magnetic moments are large 
and have opposite signs. 

For example, a valence proton polarizes, by the spin- 
spin interaction, core neutrons and these core neutrons 
give a negative contribution to the nuclear magnetic mo- 
ment (polarization of the core protons by the valence 
proton is not important). We can estimate this offset by 
considering contributions to \x from the valence and core 
nucleons. This means we have both neutron and proton 
spin contributions to (i: 



H = 9n (s z J + 9 P (s Zp ) + (l Zp ) 



(14) 



We neglected here a small contribution of the exchange 
currents into the magnetic moment. 

We want to evaluate the corrections to the valence 
model results using the very accurate experimental val- 
ues of nuclear magnetic moments. Since there are three 
unknown parameters ((s z „), (s z ), {l z )) and only one ex- 
perimental value (the total magnetic moment) available, 
to perform an estimate we need to make approximations. 

However, as we will show below, the result is not sen- 
sitive to particular approximations if we can reproduce 
the experimental magnetic moment exactly. Indeed, us- 
ing the equations above, we can present the variation of 
the magnetic moment in the following form 



5(i-5(i° = lOA5({s Zn )-(s Zn )°)-Q.56((s Zp )-{s Zp )°)} 



Sm q 
m q 
(15) 

Here 5(i° is the valence model value. For brevity we here 



Srn Q 



Sm s 



(the coefficient before 



5m s 



assume that 

in 

small anyway)). 

Let us start with the case of a valence proton. The sim- 
plest assumption is that the spin-spin interaction trans- 
fers part of the proton spin to the core neutron spin, i.e. 

((O - (0°) = -«»*,> 7 KT) and <^> = (l Zp )°. 
Then we can solve the equation for the magnetic moment 
and obtain the deviation from the valence model value 



5(i-5(i° = -0.11{fi - fi°) 



Sm a 



(16) 



To test the stability of the result, we can try different 
"extreme" assumptions. For example, if the angular mo- 
mentum exchange occurs exclusively between the pro- 
ton spin and proton orbital angular momentum, then 

this case 



5(1 - 5(i° = -0.12(/i- (i°) 



5m 



(17) 



Finally, we can try an unreasonable assumption that the 
exchange happens between the proton spin and neutron 
orbital angular momentum: (l z ) — (l z ) = —((s z ) — 

W), <0 = {**>)", (h P ) = Then 



5(i - 5(i° = -0.10(/i-M°)- 



(18) 
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We see that the results are very stable, the difference 
in the correction to the valence model is about 10 % 
only. The results for a valence neutron are similar. The 
coefficients are -0.11 for ((s Zn ) — (s Zn )°) = —((s z ) — 
(s Zp )°), -0.12 for (l Zn ) - (l z J° = -((O - (s Zn )°), and 
-0.09 for (l Zp ) - (l Zp )° = -((s z J - (s Zn )°)- 

To present the final results, we will use a physical ap- 
proximation which gives results somewhere in between 
the "extreme" assumptions. The nuclear magnetic mo- 
ments are reproduced with a reasonable accuracy by the 
RPA calculations. In the RPA approximation there are 
two separate conservation laws for the total proton j v 
and total neutron j n angular momenta (see e.g. |19j). 
We also assume that total orbital angular momentum 
(l Zn ) + (l Zp ) and total spin (s Zn ) + (s Zp ) are conserved 
(this assumption corresponds to neglect of the spin-orbit 
interaction). We repeat again that we only need these 
approximations to obtain specific numbers which are in 
between "extreme" model values. Then we can write 



{S Z }° = (S Z 
(jz p ) = {L 



,> + <»*»> 



(19) 
(20) 



where (j Zp } = / for a valence proton and (j Zp } = for a 
valence neutron. Using equations 120|) and l|19|) to elimi- 
nate (l z ) and (s Zp ) in equation l|14f> we get 



f 1 = 9n(s Zn ) + (g P - l)(s 2p ) + (j Zp ) 



M- 0%) ~ (9p~ 1 )( s z )° 



9p 



(s Zp ) = (s z )° - (s z J 



(21) 



(22) 



(23) 



We thus have taken into account both the proton and 
neutron contributions to the nuclear magnetic moment 
and can more accurately estimate how a variation in 
quark mass relates to a variation in fj,. From equation 
(|14fl we see immediately that 



<fy« = Sg n (s Zn ) + Sg p (s Zp ) 



and thus 



where 



K„ 



5fx 

9n(s Zn ) 



&9r, 
9n 



5g P 

9 P 



A' 



(24) 



(25) 



and K„ 



9 P (sz p ) 



From the definition of the g-factor for free protons and 
free neutrons, we know Sg n /g n = 5(i n / ' \x n and Sg p /g p = 
S/ipf/ip. We can now use equations (0-0 to explicitly 
relate the variation in /i to the variation in quark masses. 
Thus 



Sm„ 



(26) 



where clearly 



n q = -0.118A„-0.087A P 



k s = 0.0013A„ - 0.013A„ 



We are now in a position to evaluate the coefficients 
in specific cases. For 133 Cs, r = 7/2 + and ,u = 2.5820 
and it has a valence proton. Thus (j z ) = i = I — h 



and {s z )° = —jg- Therefore equations 
immediately give us 

{s z J = -0.103 



and 



(Sz p ) 



-0.286 



and thus 



giving 



0.152- 



0.619- 



9n 



9 P 



5u Sm„ 5m„ 

_L = 0.0358 q - +0.00824 -. 

(i m q m s 

The dependence on the strange quark mass is relatively 
weak and it is convenient to assume that the relative vari- 
ation of the strange quark mass is the same as the relative 
variation in the light quark masses (this assumption is 
motivated by the Higgs mechanism of mass generation) . 
We restore the explicit notation for the strong-coupling 
constant and conclude 



0.0441 



5{m q /A QC D) 
(m q /k QCD ) 



for 



l.TS 



Cs. 



(27) 



Values for 87 Rb, 199 Hg, 171 Yb, m Cd and 129 Xe can 
be similarly calculated and are presented in Table ITT1 

For 2 H and 3 He, the magnetic moments are pretty close 
to naive values, therefore we have not tried to improve 
the results. 

These results can be summarized in Tabic [H] A com- 
parison with the earlier results (using the valence model 
method) for the total variational relation {n q + k s ) is 
shown later in Table ITTT1 We see that 5/i/n in the nuclei 
with a valence proton is very sensitive to the core polar- 
ization effects. However, there is no such sensitivity in 
the nuclei with a valence neutron. To explain this conclu- 
sion one should note that a neutron does not give an or- 
bital contribution to the nuclear magnetic moment. The 
orbital contribution of core protons is relatively small. 
As a result, the core polarization effect changes Sfi and fi 
in a similar way, i.e. it practically does not change their 
ratio. In the nuclei with a valence neutron 



— = -0 117 ^^ - L25 ^ 6m i 
V {SzJ - 1.20(s 2p ) m q 



(28) 
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TABLE II: Variation in n incorporating the effect of non- valence nucleons in various atoms. (s z }° is the spin expectation value, 
(s z „) and {s Zp ) are (either, depending on nucleus) the valence and non- valance nucleon contributions to this spin, K„ and K v 
are defined in equation 1251 and K q and defined in equation 1261 



Atom 


X H 


2 H 


3 He 


87 Rb 


ln Cd 


129 Xe 


133 Cs 


171 Yb 


199 Hg 


(Sz)° 





1 


0.5 


0.5 


0.5 


0.5 


-0.389 


-0.167 


-0.167 


(*.»> 





0.5 


0.5 


0.124 


0.343 


0.365 


-0.103 


-0.150 


-0.151 


( S *p) 


0.5 


0.5 





0.376 


0.157 


0.135 


-0.286 


-0.017 


-0.016 


K n 








-0.172 


2.21 


1.80 


0.152 


1.16 


1.14 


K p 








0.764 


-1.47 


-0.969 


-0.619 


-0.194 


-0.173 


HLq 


-0.087 


-0.020 


-0.118 


-0.046 


-0.133 


-0.128 


0.036 


-0.120 


-0.120 


Kjs 


-0.013 


-0.044 


0.0013 


-0.010 


0.022 


0.015 


0.008 


0.004 


0.004 


Kq -\- K s 


-0.100 


-0.064 


-0.117 


-0.056 


-0.111 


-0.113 


0.044 


-0.116 


-0.116 



C. Effect of variation of the the spin-spin 
interaction 

In previous subsection we have not taken into account 
that (s Zp ) and (s Zn ) may depend on the quark mass. 
However, this dependence appears since the spin-spin in- 
teraction depends on the quark masses. Below we want 
to perform a rough estimate of this effect using the one- 
boson-exchange model of the strong interaction. 

Consider, for example, a nucleus with a valence proton. 
As above we assume that the total spin of nucleons is 
conserved and so 

<0° = (**p> + <»0 (29) 

where (s z )° is again the valence model value. It is con- 
venient to use the following notations: 

(O = b(s*)° (30) 
(s Zp ) = (l-b)( Sz )° (31) 

where b is a coefficient determined by the spin-spin in- 
teraction. We need to determine the dependence of b 
on m q / Aqcd to complete our calculations. It can be 
estimated using perturbation theory as 

where 

V aa = Vflri -ra|)Si ■ S 2 

and E Q — Ek is the spin-orbit splitting (see e.g.jl^)- As 
Trig — > 0, E — Ek remains finite and so it can only have a 
weak dependence on m q / Aqqd- The major dependence 
comes from the 7r-meson mass which vanishes in the chi- 
ral limit m q — > 0. According to review |20j the n- meson 
exchange gives about 1/3 of the spin-spin interaction. 
The other most significant contribution is given by the 



p-meson exchange. We neglect other meson contributions 
and assume that the remaining 2/3 of the spin-spin in- 
teraction is given by the p-meson. The result is not very 
sensitive to this assumption since the p— meson and other 
vector mesons have approximately the same and rather 
weak sensitivity to a variation in m q . According to Ref. 

M 

= 0>021 *n, 
m p m q 

whereas for the 7r— meson 



~ \fm q K Q cD 

so we have 

dm^ 1 Sm q 

TO?,- 2 m q 

The dominating contribution is therefore given by tt- 
meson exchange. The exchange contribution of is 
small due to the small overlap between ipp( r ) an d V'n( r ) 
and is not important. The main contribution is when a 
neutron is excited through ir° exchange into a spin-orbit 
doublet, j = Z + i to j = l — The strongest dependence 
thus originates from the 7r° pion mass. 

The momentum space representation of the nucleon- 
nucleon interaction due to a 7r-meson is 

V«(q) = gl(n ■ t 2 )(cti • q)(CT 2 ■ q) 2 1 2 

where q is the momentum transfer q = p! — p 2 , r is the 
isotopic spin and a is the Pauli spin matrix. 
We separate this into tensor and scalar parts 

((fi • q)(a 2 • q) ={{o\- q)((f 2 • q) - \a x ■ a 2 q 2 } 

+ {^l-cr 2 q 2 } (33) 
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The scalar part of the interaction we are interested in 
becomes 



scalar 



(<?) 



(n • f 2 )(o ; i • <7 2 ) - 



mi 



3 \ mj T + q z 



T 

,2 



(34) 



Fourier transformation of this and letting r 12 = |ra— ri | 
will give us the coordinate space representation 



ys 



9l 



\r l2 ) = ■^(r 1 -T 2 )(a 1 -a 2 )x 



-4nd(r 12 )+m 2 - 



(35) 



As m q — > 0, ^71- is finite and so we can neglect its depen- 
dence on m q / Aqcd- Now the strong force short range 
repulsion (implying proton and neutron hard cores) needs 
to be taken into account. It means that r 12 ^ and hence 
S(r 12 ) — 0. Nucleon core repulsion is incorporated into 
the interaction using the factor f{r 12 ) which is presented 
e.g. in Ref. H3 



and 



/>oo 

S 1 =4nm 7r / e- m ^[f(r 12 )] 2 [r 12 ] 2 dr 12 
Jo 



so that 



Thus 



dS Si 



6S 5m v Si 
S m n S 

From these definitions, we have 



and so 



SB ^Sm v 
B rn-K 



B = mi ■ S 



SS ^ 
S m-n 



Sm^ Si 

TOjr S 



Recalling that the it— meson contributes only 1/3 to the 
spin-spin interaction, we have: 



V(r 12 ) = [/(r 12 )]V(r 13 ) 



where 



/(r la ) = l-e- 1 - 1 " (1-0.68 p 2 ) , p = 

Clearly, this factor restricts nucleon interaction at very 
short ranges, with f(r 12 — 0) = 0, yet its effect is minimal 
at larger distances since f(r 12 ) « 1 for r 12 > 1 fm. It is 
this factor which results in a non-zero dependence of b 
on m q /A QC D- 

Therefore we have the effective spin-spin interaction 



V ss — Vo(r 12 )(Si ■ S2) x constant 



where 



V (r 12 ) = ml e" 



(36) 



(37) 



We can obtain the short range limit of this effective in- 
teraction (which takes into account the finite size of the 
nucleons, that is the short-range repulsion). For clarity 
in our equations, we define B, S and S±. First let 



V {r 12 )*B5(r 12 ) 



so that we have 



B = B J5(r 12 )d 3 r 12 = J V (r 12 )d\ 2 



Also let 



S 



Vo{r 12 )d 3 r 12 



An I e- m ^[f(r 12 )]\ 12 dr 12 (38) 



(o\Vss\k) 
E„ — Ei. 



and with 



Sm^ 1 5r 



2 m„ 



we see 



ffl = 2/ SA 6m* 
b 3 V S ) mjr 



1 / Si \ 8m n 

3 2 -i hr (39) 



The integrals for S and Si can be evaluated using m w 
m„o = 135MeV = 0.68 fm" 1 to give 



Si 

s 



2.17. 



This gives a small number for 4r — 2 = 0.17, therefore the 
result may seem to be unstable. It is useful to clarify this 
point using a simpler analytical model for the repulsive 
core with f 2 = 1 — exp (— kr) which gives 



Si _ 1-R 3 
S 1-R 2 



(40) 
(41) 



Any value of k > m n gives R < 0.5 and so gives a small 
difference for 4f — 2. Therefore, the small value of this 
difference does not indicate any strong instability. If we 
take k = 1.1 fm -1 (the same value of the core radius 
which we used in the more sophisticated model for / 
described above), we obtain 4+- = 2.2, i.e. practically the 
same result as above. Thus, the result does not have any 
strong model dependence. 

Using 4f = 2.17 we obtain the following 7r-meson con- 
tribution to the variation of b: 



-0.057 ^. 



(42) 
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Similarly, for the p-meson we obtained Sl/S — 3.77 and For 133 Cs, we use the values obtained earlier to get 



6b 
T 

The final estimate is 

5b 
~b 



-2.4^ 



-0.0A 



-0.11 



Sm„ 



(43) 



(44) 



total 



Note that if we assume that the spin-spin interaction is 
completely dominated by the 7r-meson exchange the re- 
sult (-0.17) would not be very different. 

Returning to the magnetic moment, we have for the 
case of a valence proton: 

M = 9n b (s z )° + (g p - 1)(1 - b)(s z )° + (j Zp ) 



6 t = S9 ILKn+ 6j E Sb 
A* 9n g P o 



K h 



(dn - 9p + l)(Sz n ) 



0.335 



giving 



^ = 0.009 ^ m i/ A Q CD ^ 



(m q /A QC D) 



All calculations for 139 La, 87 Rb, 199 Hg, 171 Yb, ln Cd 
and 129 Xe are similar to the method used for 133 Cs. The 
results are presented in Table ITTll which summarizes the 
three methods used. Method A was the first method 
discussed and used the theoretical nuclear magnetic mo- 
ment of just the valence nucleon. Method B included 
the contribution from non-valence nucleons to the nu- 
clear magnetic moment. Method C further included the 
effect of a variation in quark mass on the spin-spin in- 
teraction itself. It shows the significance of the spin-spin 
interaction on how \x varies with quark masses, with sign 
reversal for some nuclei. 



Kh 



(g n -9 P + i)(s Zn ) 



(45) 



Similarly, for the case of a valence neutron: 
A 4 = .9n(l - b)(s z )° + (g p - 1) b (s z ) c 



At 9n 



K K = 



9 P 



(.9p 



(46) 



The dependence ofSfi/fi on the spin-spin interaction (via 
the K b term) can now be seen to be quite significant. It 
depends on three values. The first is the common factor 
(dp ~ 9n — 1) = 8.41, which is large. Next, it depends on 
the value of the effective spin of the non-valence nucleons, 
which indicates the extent of the spin-spin interaction. 
It is most significant when the experimental value for /i 
deviates greatly from the valence model value. Third, 
it depends on the nuclear magnetic moment and so is 
further enhanced when dealing with nuclei with small 
magnetic moments (e.g. 111 Cd has n — —0.5949 whereas 
133 Cs has fj, = 2.582). 

We now have a modified version of equation (12(it which 
includes the term — 0.11K o to account for the variation 
of the spin-spin interaction itself: 



6 ft 
A« 



S(m q /A Q cp) 
" (mq/AQCD) 



K = -0.12K n - O.lOKp - 0.11Kb 



(47) 



Limits on the variation of the fine structure constant 
q and (rriql Kqcd) using recent atomic clock 
experiments 

We can now estimate the time dependence of the ra- 
tio of the hyperfine transition frequencies to variations in 
m q l Aqcd- The results for each atom M can be presented 
using the parameter V as defined earlier with equation 
1|11JI , with the values of K re \ and n for the atoms consid- 
ered here summarized in Table Uvl 



SV(M) _ 5{A/E) 
V(M) ~ A/E 



a 



2+K r 



A 



QCD 



me 



For two atoms, Mi and M2, the dependence of the 
ratio of the frequencies A(Mi)JA(M-2) can be presented 
as the ratio X 



X(M 1 /M 2 ) 



V(Mx) 
V(M 2 ) 



A, 



QCD 



For A( 87 Rb)/A( 133 Cs), we have 



A(Rb/Cs) = a 



0.49 / " l q 



A 



QCD 



-0.025 



(48) 



(49) 



and the result of measurements by [8j can be presented 
as a limit on the variation of X: 



1 dA(Rb/Cs) 



X(Rb/Cs) dt 
For A( 133 Cs)/A( 1 H), we have 

A(Cs/H) = a 



(-0.5 ±5.3) x 10- lb /yr 



A 



QCD 



0.109 



(50) 
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TABLE III: Comparison of results for k (see equation 1471 1 for the three methods used in various nuclei 



Atom 


87 Rb 


xll Cd 


129 Xe 


133 Cs 


139 La 


171 Yb 


199 Hg 


Method A* 


-0.074 


-0.117 


-0.117 


0.127 


0.127 


-0.117 


-0.117 


Method B* 


-0.056 


-0.111 


-0.113 


0.044 


0.032 


-0.116 


-0.116 


Method C* 


-0.016 


0.125 


0.042 


0.009 


-0.008 


-0.085 


-0.088 



* Method A considered only the valence nucleon, Method B includes the non-valence nucleons, 
Method C further includes the effect of quark mass on the spin-spin interaction. 



TABLE IV: Summary of final results showing the relative sensitivity of the hyperfine relativistic factor to a variation in 
a (parameter K re i) and the relative sensitivity of the nuclear magnetic moment to a variation in the quark mass/strong 
interaction scale m q / Aqcd (parameter ft). These values can be used in equation 11481 . 



Atom 


iH 


?H 


!He 


f 7 Rb 




IfXe 


133 p„ 
55 


139t „ 
57 ^ a 


^Yb 


80 H g 


K re l 











0.34 


0.6 


0.8 


0.83 


0.9 


1.5 


2.28 


K 


-0.100 


-0.064 


-0.117 


-0.016 


0.125 


0.042 


0.009 


-0.008 


-0.085 


-0.088 



and the result of the measurements in Ref. may be 
presented as 



dX(Cs/R) 



'X(Cs/H) dt 
For A( 171 Yb+)/A( 133 Cs), we have 



< 5.5 x IQ-^/yr. 



X(Yb+/Cs) = a 



0.67 



A 



QCD 



-0.093 



(51) 



(52) 



and the result of measurements by [T(| can be presented 
as a limit on the variation of X: 



1 



dX(Yb+/Cs) 



X(Yb+/Cs) dt 
For A( 199 Hg)/A( 1 H), we have 

X(Hg/H) = a 2 - 



(2.8 ±2.9) x 10" 14 /yr 



A, 



QCD 



(53) 



and the result of measurements by can be presented 
as a limit on the variation of X: 



1 



dX(Hg/H) 



X(Hg/H) dt 



< 8 x l(T 14 /yr 



The optical clock transition energy E(Hg) (A = 282nm) 
in the Hg + ion can be presented in the form: 



E(Hg) — const x ( 



m e e 



\F rel {Zu) 



and calculations by Ref [T|j gives 



6E(Hg) _ 5a 
£(Hg) 6 - a 



(54) 



(55) 



corresponding to ^(Hg Opt) = a 3 2 . Variation of the 
ratio of the hyperfine splitting A(Cs) is given by 



V(Cs) = a 2 



A, 



QCD 



0.009 



— • (56) 



The relative variation of the electron to proton mass ratio 
can be described by |l8j 



A 



X(m e /m p ) 
giving 

V(Cs) - a 2 ' 83 



QCD 



-0.037 



A 



QCD 



A, 



QCD 



-0.039 



(57) 



^qcdJ \-^qcdJ 

Variation of the ratio of the hyperfine splitting A(Cs) 
to this optical transition energy is given by: 



A(Opt) 



V(Cs) 



V(Kg Opt) 



a 



\ 0.009 / 

J^-) ( 2£ ) (58) 
^qcdJ \m : 1 



-0.039 



in, 



-) (59) 



^QCD/ \Aqcd J 

and the result of measurements by |l2j | can be presented 
as a limit on the variation of X: 



1 



dX(Opt) 



A (Opt) dt 



< 7 x l(T 15 /yr. 



For the Is — 2s transition in hydrogen the relativistic 
corrections are negligible, i.e. V(H Opt) = a°. Variation 
of the ratio of the hyperfine splitting A(Cs) to this optical 
transition energy is given by: 

0.009 



A(Opt) 



V(Cs) 



,2.83 



V(H Opt) 

— „,2-83 



A 



A 



QCD 



QCD 
0.039 



A 



QCD , 



(60) 
(61) 
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and the result of measurements by fl3l ] can be presented 
as a limit on the variation of X: 



1 dX(Cs/R Opt) 

AT(Cs/H Opt) It 



(3.2 ±6.3) x 10- ls /yr 



For the optical clock transition energy E(Yb) (A 
436nm) in the Yb + ion calculations by Rcf 23] gives 



M?(Yb) _ Sa 
"B(YbT " ~ 



(62) 



corresponding to V(Yb Opt) = a 88 . Variation of the 
ratio of the hyperfme splitting A(Cs) to this optical tran- 
sition energy is given by: 



X(Opt) 



V(Yb Opt) 



a 



A 



1.95 / " l 1 



A, 



QCD 





0.009 


m q 


\ m e 


-QCD; 


/ m p 


-0.039 






( m e \ 


1 


\ Aqcd J 



(63) 



(64) 



and the result of measurements by [T^j can be presented 
as a limit on the variation of X: 



1 dX(Opt) 
X(Opt) It 



= (1.2 ±4.4) x Kr ls /yr 



Other combinations have been suggested as possible 
areas of research. For A( 129 Xe)/A( 3 He), we have 



X(Xe/Ke) = a 



0.8 / m q 



A 



QCD 



0.159 



while for A( 1 H)/^( 2 H), we have 



X( 1 H/ 2 H) = 



A 



QCD 



-0.036 



(65) 



(66) 



One can use Table IFvl to predict which hyperfme tran- 
sitions will be most sensitive to a variation in a. The 
greatest effect will be seen for ratios between atoms with 
the greatest difference in values of K re i and k, especially 
if relation (JIJ were correct. Clearly, it would be best to 
test ratios of elements with opposite signs for k so that 
the effects are more pronounced. 

The effect of the spin-spin interaction is to reduce the 
sensitivity of Cs to a variation in quark mass and enhance 
the sensitivity of other nuclei such as Cd. Because the 
spin-spin interaction is so strong for Cd, with its mag- 
netic moment of -0.59 being quite different to the va- 
lence model value of -1.9, it may be quite sensitive to a 
variation in quark mass. In Ref. |6(, the importance of 
Cd was motivated by its small magnetic moment. This 
could enhance its sensitivity to a variation of the fun- 
damental constants. We have not actually obtained any 
enhancement, and the absolute value of \k\ in Cd is com- 
parable to the valence model value. However, due to the 
spin-spin interaction it has an opposite sign relative to 
some other nuclei with large \k\. For example, consider 



A(Cd)/j4(H) and A(Cd)/A(He), with each ratio involv- 
ing opposite signs for k : 



X(Cd/H) = a 



0.6 



A 



QCD 



X(Cd/He) = a 



0.6 



A, 



QCD 



0.24 



(67) 



(68) 



Note that if relation were correct, the variation of 
X may be dominated by m q / Aqcd'- for A(Cd) / A(H) , 
X(Cd/H) oc a 9 . 



III. CONCLUSIONS 

The results of this work are presented in the previous 
section. Table llVl provides one with the numbers needed 
for the interpretation of the measurements. Below, we 
would like to formulate a few conclusions which correct 
some misconceptions in the existing literature and may 
help to plan future experiments and calculations. 

• There is no such thing as a "model-independent inter- 
pretation of measurements" if one uses the valence model 
(Schmidt) values of the nuclear magnetic moments. The 
valence model cannot even guarantee the order of mag- 
nitude and sign of the effect. 

The situation may be improved by presenting the nu- 
clear magnetic moment as a linear combination of the 
neutron, proton and orbital magnetic moments. How- 
ever, even this method does not guarantee high accuracy 
since the expansion coefficients in this linear combination 
depend on the ratio m q / Aqcd- A consistent interpreta- 
tion of the measurements requires the calculation of the 
dependence of nuclear magnetic moment on this param- 
eter. 

• A small value for the nuclear magnetic moment 
does not guarantee an enhancement of the sensitivity to 
m q l Aqcd- However, a large deviation from the valence 
model value should increase the error in the calculated 
sensitivity. 

• The dependence on m q j Aqcd of the nuclear mag- 
netic moments of flRb, 55 3 Cs, 5 39 La and j^Xe is strongly 
suppressed by the many-body corrections - see Table 
IIIII One cannot guarantee high accuracy of calculation 
in this situation. However, it is probably not impor- 
tant since the suppression means that the contribution 
of these magnetic moments to the final effect of the vari- 
ation is small. The effect will be dominated by the vari- 
ation of a, m e /m p or another magnetic moment. For 
I^Cd (where the magnetic moment is small) the effect 
is not suppressed but it is of opposite sign to the valence 
model value. For ^Yb and go 9 Hg the deviations from 
the valence model are small. Naive calculations give also 
reasonable results for 2 H and ?,He. 
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